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Abstract 

A Lorentz violating modification of massless QED is proposed, with 
higher order space derivatives for the photon field. The fermion dynamical 
mass generation is studied with the Schwinger- Dyson approach. Perturbative 
properties of the model are calculated at one-loop and discussed at higher- 
order loops, showing the consistency of the model. We explain that there is 
no contradiction with the definition of the speed of light c, although fermions 
see an an effective light cone, with a maximum speed smaller than c. 



1 Introduction 

Recently, field theories with higher-order space derivatives have attracted at- 
tention, because of the improvement of graph convergence PQ, at the price of 
violating Lorentz symmetry at high energies. Ghosts are not introduced by 
this procedure, since the time derivative order remains minimal, such that 
no new pole appears in the propagator of particles. Also, theories based 
on anisotropic scaling of space and time (Lifshitz type theories) allow new 
renormalizable interactions, and for example a renormalizable exponential 
potential, in 3+1 dimensions, has been studied in j2]. Finally, a renormaliz- 
able Lifshitz type theory of Gravity has been proposed, which could lead to 
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Quantum Gravity [3|. 

A Lorentz- violating extension of massless Quantum Electrodynamics (QED) 
is discussed here, where higher order space derivatives are introduced for the 
photon field. This model has isotropic scaling in space time though, such that 
the higher order space derivatives involve a mass scale M, whose role will 
be discussed a bit further. Higher order space derivatives were already con- 
sidered in a modified Dirac equation [I] , where the resulting phenomenology 
concerning gamma ray bursts is also discussed. Higher order space deriva- 
tives could arise, for example, from quantum gravitational space time foam 
[5] . The mass M naively suppresses the effect of higher order derivatives in 
the IR, but as we will see, an IR signature of these higher orders remains: 
a fermion mass is generated dynamically. This mass, although proportional 
to M, is orders of magnitude smaller than M since it is suppressed by an 
exponentially small function of the fine structure constant. We will study 
here this mass generation, using the non-perturbative Schwinger- Dyson ap- 
proach, and we will discuss the properties of the result. We note here that 
another alternative to the Higgs model has been proposed in the context 
of anisotropic theories [5], and dynamical mass generation has been studied 
in [7] for a Lifshitz-type four-fermion model, and in [8] for a Lifshitz type 
Yukawa model. 

We mention here the first works on dynamical mass generation in QED [9], 
where a relation between the bare and dressed masses is derived, which in- 
volves a cut off A. In this context, it was shown that, in order to have a 
finite theory, the limit A — > oo implies that the bare mass should vanish, 
such that the dressed mass must be of dynamical origin. The Schwinger 
Dyson equation with a finite cut off was studied later [10] , where a critical 
value for the fine structure constant was found in order to have dynamical 
mass generation in the limit A — > oo. A summary of these results can be 
found in [TT] . 

Once the fermion dynamical mass generation is demonstrated, we study the 
perturbative quantization of the model, taking into account the fermion mass. 
We show that the model leads to the same renormalization flows as QED, 
when one considers the running scale M. The only difference is the effective 
light cone seen by fermions, since space a time derivatives get different quan- 
tum corrections. We find that the maximum speed for fermions is smaller 
than the speed of light, by terms of one- loop order, showing that the model 
is consistent with causality. This is not in contradiction with the definition 
speed of light, since the latter is obtained in the massless limit, which cannot 
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be taken for fermions, due to dynamical mass generation. The speed of light 
is given by the dispersion relation for the gauge field, which is not altered by 
the Lorentz- violating model proposed here. 

We note finally that more general higher order derivative extensions of QED 
have been presented in (12] and references therein. In these works, the au- 
thors consider Lorentz- violating vacuum expectation values for tensor fields, 
which allow the introduction of higher order derivatives for the photon field. 
They explain that the Lorentz-violating Lagrangians can be written as the 
Lagrangian for QED in a medium, and they study for example the corre- 
sponding birefringence effects of the vacuum. Our present study corresponds 
to a specific case of the latter models, where we develop quantum properties 
of a given Lorentz-violating extension of QED. 

In Section 2 we define the model and discuss some of its classical features. 
The fermion dynamical mass is then derived in Section 3, and the properties 
of the quantum theory are developed in Section 4, where we first focus on 
one-loop, and then develop general arguments for higher-order loops. The 
Conclusion presents future extensions involving non-Abelian gauge fields, and 
an alternative to the Higgs mechanism. Finally, the Appendix contains the 
details of one-loop calculations. 

2 Model 

The Lorentz-violating Lagrangian considered here is 

^ = -\ F " V (l " ^) ^ - \W (l - d » A " + ® W> (!) 

where = d^ + ieA^, and A = — did 1 = d-d (the metric used is (1,-1,- 
1,-1)), which recovers QED in a covariant gauge if M — > oo. The Lorentz- 
violating terms have two roles: introduce a mass scale, necessary to generate 
a fermion mass, and lead to finite gap equation, as will be seen further. We 
stress here that M regularizes only loops with an internal photon line, and 
that another regularization is necessary to deal with fermion loops. Also, the 
Lorentz violating modifications proposed in the Lagrangian do not alter 
the photon dispersion relation, which remains relativistic. 
No higher order space derivatives are introduced for the fermions, for the 
following reason: in order to respect gauge invariance, such terms would 
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need to be of the form 

j^z^iiD ■ 7)^ n > 2, (2) 
such that new and non-renormalizable couplings would be introduced. 



Equations of motion 

The classical equation of motion for the free gauge field is 

(l - (Drj^ - d,d u ) A» = , (3) 

and has two solutions. One is the usual plane wave, with the usual dispersion 
relation, and the other is solution of 

If the solution depends on time and one spatial direction x only, it reads 

A, = K,(t)e- M * , (5) 

where K^{t) is an homogeneous vector. If we consider a spherically symmetric 
solution, depending on time and the radial spatial coordinate r only, we find 

-Mr 

A, = JT„(t) — . (6) 

In both cases, the only regular solution occurs when = 0. As a conse- 
quence, the only physical solution of the classical equation of motion for the 
free gauge field is the usual plane wave. 



Lorentz transformation 

Finally, one can ask what the effect of a Lorentz transformation is on higher 
order space derivatives, as in the model ([T]). For this, we consider two coor- 
dinate systems (t,r) and (t',f '), with the following transformation law 

t = ^t' + v-f') (7) 
f = f± ' + 7(f]| ' + vt') , 
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where the usual notations are used. In the Lorentz transformation ([7]), we 
have 

- d 2 + A = -d 2 t , + A' with d t > = jd t + ■ V , (8) 

such that 

A' = A + ( 7 2 - 1)<9 2 + 2 7 2 {7- V<9 4 + 7 2 (tJ- V) 2 . (9) 
Using eq.([H]), the equation of motion (M 2 — A')<p(t', f ') = for a scalar field 

0(t', r ') = r) = 0o exp(iu;t — ik ■ r) , (10) 

leads then to the dispersion relation 

, v-k (v-k) 2 M 2 + k 2 
u 2 + 2u— + + = , 11 

such that 

w = ± — \/m 2 + k 2 - (v- kf/v 2 . (12) 

As a consequence, the corresponding solution is not observable, since it de- 
cays in a time of the order of the Plank time (we do not take into account 
the non-physical solution increasing exponentially in time). Therefore the 
higher order time derivatives introduced by a Lorentz transformation do not 
affect the model presented here. 



3 Gap equation 

We review here the results of [13J. The Schwinger-Dyson equation for the 
fermion propagator is [11] : 

G- 1 - G^ e = I D^WGr, (13) 

where T u , G and D^ v are respectively the dressed vertex, the dressed fermion 
propagator and the dressed photon propagator. This equation gives an exact 
self consistent relation between dressed n-point functions, and thus is not 
perturbative. As a consequence, no redefinition of bare parameters can be 
done in order to absorb would-be divergences, and for this reason one needs 
this equation to be regularized by M. 

In the studies of dynamical mass generation in QED in the presence of an 
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external magnetic field B ("magnetic catalysis" [II]), the gap equation, in the 
Lowest Landau Level approximation, is finite because of the mass scale y\eB\ 
- where e is the electric charge - which plays the role of a gauge invariant cut 
off. Another example of dynamical mass generation, however of a different 
nature, is the Debye screening for the photon, at finite temperature, which 
was found to be enhanced by a strong magnetic field [15] • 
In order to solve the Schwinger- Dyson equation f[T3j) . we consider the so- 
called ladder (or rainbow) approximation, consisting in taking the bare vertex 
instead of T u . It is known that this approximation is not gauge invariant [XT] , 
but, as we will see, the gauge coupling dependence of the dynamical mass is 
not affected by the choice of the gauge parameter £. Then we will neglect 
corrections to the photon propagator, which is given by 

D bare (cj v) = - I ^ + C — ^ (U) 

K ,P) 1+p^/M 2 \-oo 2 +p 2 V 2 -P 2 ) 2 /' 

where ( = !/£, — l,po = u and p 2 = p- p. Also, we neglect the fermion wave 
function renormalization: only the fermion dynamical mass will be taken 
into account as a correction, such that the dressed fermion propagator will 
be taken as 

u 2 - p 2 - m z dyn 

where is the fermion dynamical mass. With these approximations, the 
Schwinger-Dyson equation f[T51) - involving a convergent integral - leads to 



a 



OO fOO 



r_ I a I p dp m dyn {A + C) nfi x 

m dyn ^ j ^ au^ 1 + p2/M2 (w2 + p2) {uj2 + f + ^ , ( i > , 

where a = e 2 /47r is the fine structure constant. This equation has the obvious 
solution rridyn = 0, and potentially a second solution, which must satisfy the 
following gap equation, obtained after integration over the frequency u, 



(4 + C)a J 1 + /i 2 x 2 V VT 



x 



where \x = md yn /M is the dimensionless dynamical mass, expected to be very 
small fi << 1, and x = p/m^n. Both terms in the last equation, if taken 
separately, lead to diverging integrals, with divergences which cancel each 
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other. An integration by parts for the second term leads then to 



2tt J_ f°° ln(l + fj. 2 x 2 ) 



(4 + C)« /i 2 7 (1 + x 2 ) 3 / 2 



V V27 v^T? 
= ln(£)-^ + 0(// 2 ln//) (18) 

and the fermion dynamical mass is finally given by 

m^*Mew(-j^^y (19) 

Note that the expression (fl9|) for m^ n is not analytic in a, such that a 
perturbative expansion cannot lead to such a result, which justifies the use 
of a non-perturbative approach. A perturbative expansion would lead to the 
solution rridyn = only. 

Among the two solutions rridyn = and md yn 7^ 0, the physical system 
chooses the non-vanishing dynamical mass, in order to avoid IR instabilities, 
not favorable energetically, which would otherwise occur in the theory. 



Gauge dependence 

There is an obvious dependence on the gauge parameter (, which has a con- 
sequence on the value of m<f yn , but the important point is the non- analytic 
a-dependence of the dynamical mass, which is not affected by the choice of 
gauge: the resulting dynamical mass is of the form Mexp(— c/a), where c 
is a constant of order 1. This feature is known in the studies of dynamical 
mass generation in QED in the presence of an external magnetic field [H]. 
It has been argued though, using the pinch technique [T6j . that for the calcu- 
lation of physical quantities, which must be gauge invariant, one should use 
the Feynman gauge ( = to calculate the fermion self energy. The reason 
for this is the cancellation of terms arising from longitudinal contributions in 
the photon propagator. Nevertheless, if M is the Plank mass, the Feynman 
gauge for the dynamical mass ffl9l) leads to a too small result for the electron 
mass. But a multibrane scenario has been proposed in [17] , where a Randal- 
Sundrum warp factor enhances the dynamical mass, which then acquires the 
right order of magnitude for the electron, when M is the Plank mass. 
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4 Perturbative analysis of the model 



Taking into account the fermion mass md yn , we discuss here the perturbative 
quantization of the model, first at a one-loop, and we give then the general 
arguments at higher order loops. 

4.1 One-loop properties 

Fermion kinetic term 

We calculate in the Appendix the one-loop quantum corrections for the 
fermion kinetic term, which are different for time and space derivatives. We 
find the corrections 

# ((1 + Z )d i° - (1 + Z x )d ■ 7) V , (20) 

with 

z " = "^('"(m) +41n2 - 2 ) +o(p 2 Hi/v)) 

Zl = -s( ln Q) + y-f h2 ) + (" 2 ™- < 21 > 

Note that the dominant term, proportional to ln(l//i), is the same for Z 
and Zi, since in the Lorentz symmetric situation (M — > 00 for fixed fermion 
mass), we have Z = Z\. Also, the coefficient — a/(2n) in front of the domi- 
nant term ln(l///) is the coefficient found in QED in 4 — e dimensions 

Z QED = Z QED = fi (22) 

° 1 27T6 K J 

An important remark should be made here: because of the result ffTOl) . the 
ratio fi is actually finite in the limit where M — > 00, and one could think 
that no counter term is necessary to absorb terms proportional to ln(l/yu). 
But it is in fact necessary, in order to respect the loop structure. Indeed, if 
one keeps ln(l///)-terms in the renormalized theory, the would-be one-loop 
correction would become a tree-level one: 




Therefore, provided one treats ln(l//i) as a 1/e term in dimensional regular- 
ization, one gets the usual QED one-loop structure. 
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Maximum speed for fermions 

After redefinition of the bare parameters in the minimal substraction scheme, 
where only the term proportional to ln(l//x) is absorbed, the fermion disper- 
sion relation is 

(l - ^[2 In 2 - 1])V = (l - ^[25/9 - (10/3) ln2]) V + m% n , (24) 

and the product of the fermion phase velocity v$ and group velocity v g is 
then 

2 <jJ d<jj . . 

v = v+v g = -- (25) 

1 - (a/vr)[25/9 - (10/3) ln2] x 

1 - (a/7r)[21n2 - 1] 

2a /34 16 , \ , nl 9s 
= l- — l K j- Y \n2j+0(a 2 ) <1, (26) 

which shows that the effective light cone seen by fermions is consistent with 
causality. If we take a ~ 1/137, we obtain for the fermion maximum speed 

v ~ 1 - 1.9 x 10" 4 . (27) 



Gauge invariance and speed of light 

We calculate in the Appendix the one-loop vertex for vanishing incoming 
momentum, and we find the corrections 

= Z ie f, (28) 

such that gauge invariance is respected, since quantum corrections to the 
vertex are identical with corrections to the fermion kinetic term, for time 
and space components independently. As a consequence, as in usual QED, 
quantum corrections to the coupling constant are given by the polarization 
tensor, and arise only from the wave function renormalization of the gauge 
field. Since the one-loop polarization tensor doesn't include any internal 
photon line, the one-loop running of the coupling constant is therefore also 
the same as in QED. 

Because Z ^ Z x and therefore the effective light cone seen by fermions 
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involves the speed v < 1, one might see a problem with the definition of 
speed of light. We argue here that it is not the case, because of dynamical 
mass generation for the fermion. Indeed, the speed of light c is defined by 

c = lim ^- , (29) 

for finite momentum p and frequency u. But because the fermion is always 
massive, m = m dyn ^ 0, the limit f )29|) cannot be taken, and the result f )27|) 
is not in contradiction with the speed of light. Such a conclusion was already 
obtained in [8] for a Lifshitz-type Yukawa interaction. 

The speed of light is given by the gauge field dispersion relation, which is not 
modified at one-loop. We show in the next subsection that this conclusion 
still holds at higher order loops. 



4.2 Higher order loops 

Perturbative expansion 

We have seen at one- loop that, in order to recover the QED renormalization 
structure, one needs to absorb terms proportional to ln(l//i) in the redefini- 
tion of bare parameters. If we assume that, after including the appropriate 
counter terms, we have a finite theory at n — 1 loops, a graph calculated at 
n loops will be of the form 



G (n) = p(n)Ml + Q(n) ln (T) + finite 



(30) 



where, compared to (n — 1)— loop graphs: 

pi n ) contains graphs with an additional fermion loop, which needs to be 
regularized dimensionally; 

Q( n ) contains graphs with an additional photon line, which are regularized 
by M. 

The n-loop counter term then needs to absorb terms proportional to both 
1/e and ln(l//i). The contribution to the renormalization group flow is then 

M ~dM I = ^ + Q{n) + °^ ln ^ ' (31) 
and corresponds to the usual result obtained in QED, since we have 



G [ qL = (P (n) + Q (n) ) — + finite , (32) 



M' 



r QED f 
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where the finite terms differ from those in eq. (l30p . We stress here that the 
partial derivative with respect to M should be taken at fixed fermion mass. 

Speed of light 

If one considers two-loop properties of the model or higher orders, the po- 
larization tensor is affected by Lorentz violation, unlike in the one-loop case, 
and it is necessary to make sure that the model remains consistent, especially 
as far as gauge invariance and speed of light are concerned. 
From two loops and above, the field strength gets different corrections for 
time and space derivatives, and we obtain 

2(1 + Y )F 0i F oi + (1 + V ; I/-;,/-"-' , (33) 

where Yq and Y\ represent the finite quantum corrections to the operators 
and Fj* respectively, after absorbing the regularization terms proportional to 
1/e or ln(l//i). In order to obtain corrections proportional to the Lorentz 
scalar F^F^, which is necessary to recover gauge invariance and the speed 
of light c = 1, we rescale the time coordinate and the component A as: 



t — > — and An — > kAcx where k = \ / . (34) 

One can easily see then, that this rescaling is consistent with gauge invariance 
of the fermion sector: 

i) {id - eA ) -y°ip « ijj {id - eA ) j°ip (35) 

The factor k, which does not appear in the space components of the covariant 
derivatives, will then contribute to the maximum speed for fermions, together 
with the corrections to the fermion kinetic terms, as explained for the one- 
loop case. A final identical rescaling for all the gauge field components, by 
the factor a/1 + Y\, will lead to the redefinition of the coupling constant. 



5 Conclusion 

We have shown that Lorentz- violating higher-order derivatives for an Abelian 
gauge field lead to a fermion dynamical mass, independently of the strength 
of the coupling constant, and that it also leads to a consistent quantum the- 
ory. The order of magnitude of the fermion dynamical mass is too small for 
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the electron, though, if one identifies M with the Plank mass, and this is a 
problem which could be solved with the introduction of a non-Abelian gauge 
field, as explained below. 

The main difference with QED is the maximum speed for fermions, which is 
smaller than the speed of light. This implies that the fermions see a different 
effective light cone, but it is not in contradiction with the definition of speed 
of light, since the fermion is always massive. 

The next step is to extend this study to non-Abelian gauge theories. It is 
known that vector fields coupled to the axial current acquire a mass when 
a fermion condensate forms [IB] , as the result of the massless bound state 
excitation, which implies a pole in the polarization tensor for vanishing ex- 
ternal momentum. The residue of this pole is the vector mass squared. In 
the framework of our present model, the fermion condensate can occur, due 
to Lorentz-violation in the Abelian sector, and then generate dynamically 
the non-Abelian vector masses. We note that, in order not to break the 
non-Abelian gauge symmetry at the classical level, a Lorentz-violating non- 
Abelian dynamics would involve higher order covariant-spa.ce derivatives, and 
would therefore introduce new interactions, which are not renormalizable. 
For this reason, we are planning to leave the non-Abelian sector unchanged, 
and introduce Lorentz-violation in the Abelian sector only. 
Another consequence of the coupling of fermions with several gauge fields will 
be to change the order of magnitude of the fermion dynamical mass. Indeed, 
as was already seen in [TjJ], where a similar study was done, the contributions 
of different gauge fields lead to a partial cancellation in the exponent in the 
exponential of eq. f|T9|) . giving a phenomenologically more realistic dynamical 
mass. 

Acknowledgments This work is partly supported by the Royal Society, 
UK. 

Appendix: One-loop Feynman graphs 

We calculate here one-loop corrections to the model ([T]), assuming that the 
fermion has mass m, which is dynamically generated. The metric we use 
throughout is (+, — , — , — ). 
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Fermion kinetic terms 

The fermion wave function renormalization is obtained by differentiating the 
fermion self energy with respect to the external momentum u, k, and then 
set v — 0, k — 0. The fermion self energy is 

£ (1) (z/,£) (36) 
f duj f d 3 p 7 M ( ^+ ^ + m) 7 ^ 



2vr J (2tt) 3 (l + p 2/M 2 )(w 2 - p 2 )[(w + v) 2 - (p + £) 2 - m 2 ] 

_ '(27T 4 ) J ^ 7 1 + P 2 /M 2 ( W 2 _ ^2) [(^ + z/ )2 _ (£ + fc)2 _ m 2] • 

We write 

E« (i/, £) = Z z/ 7 ° - Ziife • 7 + C(A; 2 ) , (37) 

— * 

and a derivative with respect to /c p gives, after setting v = 0, fc = 0, 
lie 2 f d 3 p f duo ( 7 P 2p p 



(27r) 4 J l + p 2 /M 2 J u 2 — p 2 \u 2 — p 2 — m 2 (u 2 — p 2 — m 2 ) 2 
For p = 0, we obtain, after a Wick rotation, 

p 2 dp f dou —oo 2 + p 2 + 



e 



2tt 3 y 1 + p 2 /M 2 y w 2 + p 2 (w 2 + p 2 + m 2 ) 2 
and integration over the frequency u gives 



(39) 



Z n = 



e 2 [°° x 2 dx ( 1 , A 

J 1 + /i 2 X 2 V ^ / 



2tt 2 y l + /i 2 x : 
2a /l _ 1 



= -^( i ln(l/ / i) + ln2--)+0( / i 2 ln(l//i)), (40) 



where \x = m/M << 1 and x = p/m. 

For p = i, we obtain, after a Wick rotation, 



e 2 /" p 2 dp /" do; w 2 +p 2 /3 + m 2 
1 = ~ 2^3 J 1 + p 2 /M 2 J co 2 + p 2 (to 2 +p 2 + m 2 ) 2 ' ^ j 
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where the constant 1/3 in factor of p 2 in the numerator comes from the 
Galilean symmetry in the 3-dimensional space. The integration over u gives 

e 2 f°° x 2 dx ( 2x 1 2x 2 /3 - 1 x 2 /3 



2tt 2 y 1 + p 2 x 2 V 3 x v^Tl? 2 (1 + x 2 fl 2 

la fl 25 5 
— -m(l/u) H 



2rv /1 25 5 

' .]b(1/^) + — --ln2 ) +0(/i 2 ln(l//i)) . (42) 



Vertex 

The one-loop correction to the vertex is 



duj f d 3 p 1 7 M ( p l +m) , -f p ( p / +m)'~f^ 



2tt 7 (27r) 3 l+p 2 /M 2 (w 2 -p 2 )(w 2 -p 2 -m 2 ) 2 

e 3 /" p 2 dp f , 27 p (w 2 — p 2 ) — Ap p 4— 2m 2 Y . . 

duj ,..■>. „9 „ono/,o ox • (43) 



47T 3 J l+p 2 /M 2 J (uj 2 - p 2 - m 2 ) 2 (uj 2 - p 2 ) 

For p = 0, we obtain after a Wick rotation 

^ n n e 3 /" p 2 <ip /" , -w 2 + p 2 + m 2 



{1) 1 2n 3 J l+p 2 /M 2 J (u 2 + p 2 )(u 2 + p 2 + m 2 ) 2 

= Z Q erf . (44) 

For p = i, we obtain after a Wick rotation 

j o e 3 f p 2 dp f cu 2 +p 2 /3 + m 2 

{1) ~ ^ 2^ J l + p 2 /M 2 J UJ (u 2 +p 2 )(u 2 + p 2 + m 2 ) 2 

= Z x eff . (45) 
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